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Abstract 

We discuss electron acceleration and heating during collisionless magnetic reconnection by 
using the results of implicit kinetic simulations of Harris current sheets. We consider and 
compare electron dynamics in plasmas with different (3 values and perform simulations up to 
the physical mass ratio. We analyze the typical trajectory of electrons passing through the 
reconnection region, we study the electron velocity, focusing on the out-of-plane velocity, and 
we discuss the electron heating along the in-plane and out-of-plane directions. 
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I. INTRODUCTION 



Collisionless magnetic reconnection plays an important role in energetically active pro- 
cesses in plasmas 0]. Magnetic reconnection takes place in plasmas characterized by 
different values of (3. Theoretical, observational, and experimental results show that recon- 
nection is present in the geomagnetic tail 3JL where local (3 » 1; in the Earth's magne- 



topause 0|, where (3 « 1; in laboratory 
in astrophysics! plasmas, such as extragalactic jets 
Nuclei (AGN) yj, where (3 < 1. 

During magnetic reconnection, magnetic energy is converted into kinetic and thermal 
energy of electrons and ions. In fact, electron heating and acceleration are signatures of 
magnetic reconnection. In the magnetotail, bursts of energetic electrons have been at- 
, ribut ed to _ ti o„ Hflfland there have bee n recent di.ee, _ne„ tS of 

n 

electron acceleration during magnetic reconnection [lflj . Production of runaway electrons 
during sawtooth instabilities and disruptions is associated with magnetic reconnection in 
tokamaks j^j. In solar flares, x-ray observations indicate that a large fraction of the total 



10|l . in the solar corona plasma (h | , 



and 



1J] and flares in Active Galactic 



energy is released in accelerated electrons 
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231 ] . The observed synchroton radiation 



in extragalatic jets is thought to be generated by reacceleration or in-situ acceleration of 



electrons due to magnetic reconnection 
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141 ] . It has been proposed that the detection of 



hard x-ray and 7-ray from AGN is due to the presence of electrons accelerated by magnetic 
reconnection [151. 



Electron dy namics in the reconnection region 
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test particle theor 



and kinetic simulations [35 
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rave been studied using analytical argu- 
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self-consistent fluid 



The aim of the present paper is to study the electron dynamics near the reconnection 
region with self-consistent kinetic simulations of high and low (3 plasmas. The plasma (3 is 
varied from very large ((3 » 1) to small ((3 < 1) values and systems are simulated with 
an ion/electron mass ratio up to the physical value (mj/m e = 1836). We consider two- 
dimensional reconnection in Harris current sheet configurations triggered by an initial 
perturbation We introduce a guide field to reduce the plasma (3 and eliminate the 

null field region at the current sheet. To perform kinetic simulations, we use CELESTE3D 
[40! 41 . 42] ]. an implicit Particle-in-Cell (PIC) code, which models both kinetic ions and 



electrons while allowing simulations with higher mass ratios. 

We show that both the plasma j3 and the mass ratio strongly affect the electron dynamics. 
In low-/3 plasmas, the electron meandering orbits present in high (3 plasmas disappear. We 
focus on the out-of-plane electron velocity, which remains localized in the high-/? case, and 
which is sizable also far from the reconnection region in low (3 plasmas [36]. A strong 
influence of the mass ratio on the out-of-plane velocity is shown by the simulations and 
the relevant scaling law is deduced. We show that the heating process is non-isotropic in 
presence of a guide field; in particular, the particles are preferably heated in the out-of- 
plane direction. This anisotropy contributes to the break-up of the frozen-in condition for 
electrons and allows reconnection to happen. 

The paper is organized as follows. Section II describes the physical system and the 
simulations. Section III presents the results of the simulations, focusing on the typical 
trajectory of the electrons that pass through the reconnection region and the evolution of 
the electron fluid velocity and temperature during reconnection. 



II. PHYSICAL SYSTEM 

We consider a two-dimensional Harris current sheet in the (x, z) plane [3^, with an initial 
magnetic field given by 

BqO) = Bq tanh(z/A)e z + B y0 e y (1) 
and plasma particle distribution functions for the species s (s = e, i) by 

/„,(*, v) = no sech 2 (,/A) (^r) -p { - [v 2 x + (v y - V s ) 2 + v 2 z ] } 
( m s \ 3/2 ( m s v 2 \ 

We use the same physical parameters as the GEM challenge 3t|. The temperature ratio 
is Ti/T e = 5, the current sheet thickness is A = 0.5<ij, the background density is = 0.2n , 
and the ion drift velocity in the y direction is V^o = 1.67V^, where Va is the Alfven velocity, 
and V e o/Vio = —T eQ /T iQ . The ion inertial length, di = c/u P i, is defined using the density no- 
We apply periodic boundary conditions in the x direction and perfect conductors in the z 
direction. 
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The standard GEM challenge parameters model reconnection in high (3 plasmas. To 
model low f3 plasmas, it is possible to consider either a entirely different equilibrium [43|, 
or one may introduce a guide field in the standard Harris sheet equilibrium. Herein, we 
follow the second approach and we introduce a guide field B y = B y0 , with a spatially 
constant value at t = 0. The simulations are performed with different mass ratios, ranging 
from mi/m e = 25 (standard GEM mass ratio) to the physical mass ratio for hydrogen, 
mi/m e = 1836. Following Birn et al. j^, the Harris equilibrium is modified by introducing 
an initial flux perturbation in the form 

A y = —A y0 cos(2nx/L x ) cos(nz/L z ) (3) 

with A y0 = 0.1B Q c/ujpi, which puts the system in the non-linear regime from the beginning 
of the simulation. 

The simulations shown in the present paper are performed using the implicit PIC code 
CELESTE3D, which solves the full set of Maxwell- Vlasov equations using the implicit mo- 



ment method 
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421 ]. The implicit method allows more rapid simulations on ion length 
and time scales than are allowed with explicit methods, yet retains the kinetic effects of both 
electrons and ions. In particular, the explicit time step and grid spacing limits are replaced 
in implicit simulations by an accuracy condition, vth, e At < Ax, whose principal effect is 
to determine how well energy is conserved. In the simulations shown below, we typically 
choose uj ce At « 0.5, Ax/di = 0.4, and Az/d{ = 0.2. 

Previous work on magnetic reconnection performed by CELESTE3D have proved that 



results from our implicit code match well the results of explicit codes |42j. Implicit simula- 
tions allow one to model physical mass ratios, with which it is possible to distinguish scaling 
laws associated with different break-up mechanisms Q]. CELESTE3D has also been em- 
ployed in a comprehensive study of the physics of fast magnetic reconnection, in plasmas 
characterized by different (3 values |45J]. 



III. RESULTS 



We have performed a set of simulations, using different mass ratios (mj/m e = 
25, 180, 1836) and introducing different guide fields: B y0 = 0, with (3 = oo at the center 
of the current sheet; B y0 = B , with (3 = 1.2; and B y0 = 5B , with (5 = 0.048. We note that 
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a guide field changes drastically the magnetic configuration of the system as the X point is 
no longer a null-point, as it is in the B v q = case. 

The dynamics of magnetic reconnection in plasmas with different values of B have been 
pointed out and summarized in a previous paper |45]]. Figure 1 shows the reconnected flux, 
A\I>, defined as the flux difference between the X and the O points for all the simulations 

n 

considered |45J . Even though an initial perturbation is applied, reconnection proceeds slowly 
during an initial transient phase (which lasts approximatively until tio ci ph 10), when the 
system adjust to the initial perturbation. Subsequently, reconnection develops rapidly until 
the saturation level is reached. Both the reconnection rate and the saturation level decrease 
when the guide field is increased. All the simulations show a similar evolution. The mech- 
anism which breaks the electron frozen-in condition is provided by the off-diagonal terms 

nn 

of the electron pressure tensor for all the guide fields considered [45, 46]. The reconnection 
rate is enhanced by the whistler dyna mics in hi gh B plasmas, and by the Kinetic Alfven 
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l|, provided that B > m e /mi. When 



Waves dynamics in low B plasmas (4 
8 < m e jm il fast reconnection is not possible 

Theoretical results and kinetic simulations |44j, |54j show that with B y0 = 0, electrons 
flow towards the X point along the z direction, where they are demagnetized in a region 
corresponding to the meandering length, d ze . There they are accelerated by the reconnection 
electric field, E y , in the y direction. The electrons are then diverted by the B z field, gaining 
an outflow velocity in the x direction, and beco ming remagnetized at the meandering length, 
d xe . The meandering lengths are defined as j^J, S| 



m P T K 



e 2 (dB z /dx) 



1/4 



d. 



m P T P , 



(dBJ 



while the maximum inflow and outflow velocities scale as 4J, 



i/i 



(4) 



Am e T e (dB z /dx) 



1/4 



4m e T e (dB x /dz) 



i/i 



(5) 



In the reconnection region, the electrons are unmagnetized and follow complex meander- 
ing orbits, which result in a non-gyrotropic electron distribution function and in off-diagonal 
terms of the electron pressure 55| . 

When a guide field is introduced, the meandering orbits disappear. Analytical estimates 
of the guide field at which this happens are given by 3, |^. Nevertheless, the diagonal 



components of the electron pres sure tensor are unequal, which contributes to the presence of 
off-diagonal pressure termsQQ. These tenn 8 stiH eo— th e br ea k - U p m ee hanism of 
the frozen-in condition 45, 46]. In the reconnection region, the electrons flow across the field 
line in the (x, z) plane, while performing a Larmor motion around the out-of-plane magnetic 
field. Far from the reconnection region, the guide field causes additional components of the 
E x B force, which modify the ion and electron motion and cause asymmetric plasma flow. 

Below, we describe in detail the typical trajectory of an electron passing through the 
reconnection region in high and low (3 plasmas. Then, we focus on the electron fluid velocity, 
in particular on the out-of-plane velocity. Finally, we consider the electron distribution 
functions to evaluate the electron temperature. 



A. Electron trajectory 

Figure 2 shows a typical trajectory of an electron passing through the reconnection region 
in the case B y0 = 0, and Fig. 3 shows the history of its velocity and kinetic energy. 

In Fig. 2, the initial position of the particle is denoted by a plus sign, and its position 
at selected time steps by circles. 'X' marks the position of the X point. Note that periodic 
boundary conditions are applied in the x direction, as the behavior of the trajectory shows 
(the particle exits from the left and re-enters from the right). At the beginning, the electron 
is tied to a magnetic field line (magnetic field lines run mainly along the x direction). Near 
the reconnection point, the electron decouples from the magnetic field and moves along the 
z direction, reaching the X point at tu C i « 15. The particle trajectory is meandering in the 
unmagnetized region. The outflow from the reconnection region takes place as soon as the 
electron reaches a region with stronger magnetic field. Then, the electron couples again to 
a magnetic line surrounding the O point, and starts again its gyration orbit around it. 

In Fig. 3, all components of the particle velocity and the kinetic energy are plotted as 
a function of time. Initially, the electron is flowing along the x direction, with a Larmor 
motion mostly in the (y, z) plane, which is responsible for the high frequency oscillations 
of the velocity (the magnetic field line is mostly directed along x). During this phase, 
the electron kinetic energy is almost conserved. When the electron decouples from the 
magnetic field line as it crosses the reconnection region, it is accelerated by the reconnection 
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electric field in the y direction. This acceleration transfers magnetic field energy to electrons 
during magnetic reconnection. As Fig. 3 shows, while the electron is unmagnetized, the 
kinetic energy of the electron increases remarkably, showing high frequency oscillations due 
to the acceleration and de-acceleration by the electric field. When the electron leaves the 
reconnection region and again couples to the magnetic field, motion in y becomes Larmor 
motion with a bigger radius, and the velocity directed along y gained in the reconnection 
region is lost. The particle couples to a magnetic field line that surrounds the O point and 
starts to flow along it. 

In Fig. 4 the electron trajectory is traced for a low (3 plasma, with a strong guide field, 
B yQ = 5B . Initially, the particle flows along the magnetic line, which is mainly directed 
along the y direction, and executes a Larmor motion mostly in the (x, z) plane. The particle 
then accelerates towards the X point, crosses magnetic lines in the (x, z) plane, and gains an 
out-of-plane velocity which increases its kinetic energy (see Fig. 5). The particle still couples 
to the magnetic field and executes a Larmor orbit around the guide field. Meandering orbits 
are not present. In contrast to the case with B y0 = 0, the electron maintains its y velocity 
even when far from the reconnection region because now the gyration is in the (x, z) plane 
around the ^/-directed guide field. Finally, the electron drifts along a magnetic field line 
around the O point maintaining a still significant y velocity which decreases slowly because 
of the interactions of the electron with the non-drifting plasma background. 

The presence of the guide field changes the nature of electron acceleration: without guide 
field, the y velocity is lost while in presence of guide field is retained even far from the 
reconnection region 36]. 



B. Electron fluid velocity 



When ByQ = 0, both kinetic simulations and theoretical results [44. 54] show that the 
electrons are demagnetized at the electron meandering distance [see Eqs. (jl}] and have an 
inflow and outflow velocity given by Eqs. (J5|) . The scaling laws of the dimensions of the 
reconnection region and of the inflow and outflow velocity, based on the electron pressure 
as a break-up mechanism as derived in Ref. 3], have been verified up to the physical mass 



ratio 



44j . In the presence of a guide field, new components of the E x B field arise, and the 
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electron in-plane motion has been described in Refs. [45 

Here, we focus on the electron out-of-plane velocity. In Fig. 6, the velocity along the axis 
z = is depicted with B y0 = for three different mass ratios, mj/m e = 25, 180, 1836. The 
maximum out-of-plane velocity increases with the mass ratio. We note that the out-of-plane 
velocity is sizeable only in the reconnection region. 

The out-of-plane velocity can be estimated for B y0 = 0. The electron lifetime in the 
reconnection region, r, from Eqs. (J3HHJ is approximately 

/m 2 T 2 \ 1/4 

r w d xe /v xe + d ze /v ze oc ^ J (6) 

As the magnetic field is negligible in the reconnection region, the electrons are freely accel- 
erated and the out-of-plane velocity can be estimated as 

Vy S3 — EyT (7) 

and, using Eq. (jBJ), it follows that 

Vy OC V t h,e (8) 

Since the temperature of the electrons is the same for all mass ratios, it follows that the 
electron out-of-plane velocity scales with 1/ ^Jm e . The results presented in Fig. 6 fit well 
this scaling law, as is shown in Table I. 

Figures 7 and 8 consider the effect of the guide field on the out-of-plane electron velocity. 
As the guide field allows the particles to flow more easily in the out-of-plane direction, 
the peak velocity increases remarkably when the guide field becomes stronger, Fig. 7. 
Moreover, the presence of the guide field changes the general pattern of the out-of-plane 
velocity, as shown in Fig. 8. When B y0 = 0, the out-of-plane velocity is sizeable only 
near the reconnection region, where the electrons are accelerated by the electric field. The 
out-of-plane velocity is lost when the electrons become again magnetized and are diverted 
by the B z field. In presence of the guide field, the electrons maintain their y velocity when 
they leave the reconnection region and orbit around the O point. Note that this conclusion 
is further supported by the analysis of particle orbits performed in the subsection above 
(Ill.b). 

As a final remark, we note that the out-of-plane velocity evolves during magnetic recon- 
nection, as is shown in Fig. 9. For all the guide fields considered, the electron velocity 
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increases while reconnection proceeds (the evolution of the reconnected flux in these simu- 
lations is presented in Fig. 1). After reconnection saturates, the reconnection electric field 
vanishes and the electrons are no longer accelerated. The out-of-plane velocity in the former 
reconnection region decreases abruptly. 



C. Electron temperature 

The evolution of the electron temperatures in the reconnection region, T xe , T ye , and T ze , 
are plotted in Fig. 10 for the three different guide field strengths. We note that T xe , T ye , 
and T ze are defined as the second moment of the distribution functions of the x, y, and z 
velocity 0. 

In the zero guide field case, the T xe and T ye evolution is similar (see Fig. 10a). T xe 
increases because the positive and negative outflow velocity causes an increases spread in 
the x velocity. The heating in the y direction is due to the electric field which, besides 
accelerating the electrons, spreads out their velocity, reflecting the variation in electron 
residence time in the reconnection region, which depends on their in-plane inflow and outflow 
velocity, and thus are accelerated by different amounts. After reconnection saturates, the 
heating process stops and electrons tend to thermalize, causing an increase in T ze . We note 
that the total energy of the system is conserved during the simulation within an error of the 



order of 4% 



45|. 



When the guide field is introduced, both T xe and T ze remain almost constant at the 
initial level during the reconnection process, while T ye increases remarkably. The guide 
field introduces an higher electron mobility in the y direction. Thus, the electron can be 
accelerated by the electric field more than in the B y0 = case along the y direction, and 
the y velocities spread out more, while T ye increases. 

The anisotropy in the electron temperature contributes to the break-up of the frozen-in 
condition. In fact, in the presence of a guide field, the difference between the diagonal terms 



of the electron pressure tensor contributes to t 



for the break-up of the frozen-in condition 



le off-diagonal terms, which are responsible 



P dv B 

1 xye — r\ D ^ We zze > v I 



u ce dz B y 
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Since P yye - P zze = n e {T ye - T ze ) and P yye - P zze = n e (T ye - T ze ), the importance of the 
anisotropy in the electron temperature is evident. 

IV. CONCLUSIONS 

In the present paper, the electron dynamics during magnetic reconnection has been stud- 
ied by showing and discussing results of kinetic simulations of Harris current sheets. Simu- 
lations with different plasma (3 and different mass ratio have been considered. 

By varying the guide field, we have been able to model reconnection in systems such 
as the magnetotail (B y0 = 0), the magnetopause (B y0 = B ), laboratory and astrophysical 
plasmas (B y0 = 5B ). 

By studying the typical electron trajectories, we have shown that, when the plasma (3 is 
decreased, the electrons mainly perform Larmor motion around the guide field even in the 
reconnection region, and that meandering orbits disappear. In all the cases, electrons are 
accelerated by the reconnection electric field along the y direction and their velocity increases 
with the guide field. Moreover, the out-of-plane velocity increases during reconnection. In 
high (3 plasmas, the out-of-plane velocity is sizable only in the electron reconnection region. 
With a guide field, the out-of-plane velocity is globally relevant. The mass ratio has a 
strong influence on the out-of-plane velocity and the scaling law of interest is derived. The 
study of the electron temperature in the reconnection region has shown a strong heating 
anisotropy in presence of a guide field, which contributes to the break-up of the electron 
frozen-in condition. 

In closing, we note that we plan to develop the present work in two directions. First, we 
plan to introduce the relativistic equations of motion in CELESTE3D, in order to represent 
better the electron physics when relativistic effects become important. Second, an experi- 
mental setup has been built at the Los Alamos National Laboratory to study reconnection 
experimentally in plasmas with different (3 We plan to compare our simulation results 
with the experiments. 
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Fig. 1 (from Ref. |45j): Reconnected flux (normalized to B c/u P i), for mj/m e = 25 
(a), mi/m e = 180 (b), and mi/m e = 1836 (c), and B y0 = (solid line), B y0 = B 
(dashed line), B y0 = 5B (dotted line). 

Fig. 2: Electron trajectory in the (x,z) plane, for B y0 = and mj/m e = 25. The 
position of the particle at different times is marked by circles, the starting position by 
a plus sign. The position of the X point is denoted by the x-mark. Note the periodic 
boundary conditions in the x direction. 

Fig. 3: Velocities v x (a), v y (b), v z (c) and kinetic energy (d), as a function of time, 
for the electron whose trajectory is represented in Fig. 2. 

Fig. 4: Electron trajectory in the (x,z) plane, for B y0 = 5B and mj/m e = 25. The 
position of the particle at different times is marked by circles, the starting position by 
a plus sign. The position of the X point is denoted by the x-mark. Note the periodic 
boundary conditions in the x direction. 

Fig. 5: Velocities v x (a), v y (b), v z (c) and kinetic energy (d), as a function of time, 
for the electron whose trajectory is represented in Fig. 4. 

Fig. 6: Electron out-of-plane velocity at z — 0, when A\l/ « 1, for the simulations with 
B y0 = and mj/m e = 25 (dashed line), mj/m e = 180 (dotted line), and mj/m e = 1836 
(solid line). 

Fig. 7: Electron out-of-plane velocity at z — 0, when A\l/ w 1, for the simulations with 
B y0 = B (a) and B y0 = 5B (b), m,i/m e = 25 (dashed line), mi/m e = 180 (dotted 
line), and mi/m e = 1836 (solid line). 

Fig. 8: Electron out-of-plane velocity when w 1, for the simulations with B y0 = B 

(a) and B y o = 5Bq (b), for the simulations with mj/m e = 180. 

Fig. 9: Evolution of the average out-of-plane electron velocity, v y in the reconnection 
region, for the simulation with mi/m e = 25, and guide field B y0 = (a), B y0 = B 

(b) , and B y0 = 5B (c). 
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• Fig. 10: Evolution of the electron thermal velocity v t h x ,e (solid), v t h y ,e (dashed), and 
Vthz,e (dotted), in the reconnection region for the simulation with mi/m e = 25 and 
B y o = (a) and B y0 = 5B (b). 

Table I. Comparison between the simulation results and the scaling law in Eq. (JHJ) 



Ratio 


Simulation 


result Scaling law 


v y (m;/m e =180) 


2.9 


2.7 


v y (mi/m e =2S) 


v y {m,i 1 'm e =1836 


7.0 


8.6 


Vy(rrii/m e =25) 
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